Gaussian transformations and distillation of entangled Gaussian states 
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We prove that it is impossible to distill more entanglement from a single copy of a two-mode 
bipartite entangled Gaussian state via LOCC Gaussian operations. More generally, we show that 
any hypothetical distillation protocol for Gaussian states involving only Gaussian operations would 



be a determ inistic protocol. Finally, we argue that the protocol considered by Eisert et al. | quant 



ph/0204052| is the optimum Gaussian distillation protocol for two copies of entangled Gaussian 
states. 

PACS numbers: 03.67.-a, 42.50.Dv 



Quantum entanglement is a key ingredient of many 
protocols for quantum information processing such as 
quantum teleportation [jjj or quantum cryptography ||. 
Usually, the entangled particles are distributed among 
two distant parties traditionally called Alice and Bob. In 
practice, the transmission channel used for this distribu- 
tion is always noisy and imperfect, which prevents Alice 
and Bob from sharing a maximally entangled state even 
if Alice can prepare such state locally in her lab. Fortu- 
nately, the errors introduced by noisy quantum channels 
can be overcome by the so-called entanglement distilla- 
tion protocols, by which Alice and Bob can extract from a 
large number of weakly entangled mixed states a smaller 
number of highly entangled almost pure states || f| . 

Recently, a great deal of attention has been devoted 
to the quantum information processing with continuous 
quantum variables and continuous-variable analogues of 
various protocols developed originally in the framework 
of discrete quantum variables have been established. Re- 
markably, linear optics, parametric amplifiers, and ho- 
modyne detectors suffice for implementation of many of 
these protocols including continuous-variable teleporta- 
tion 0, cryptography ||, and cloning Q. However, one 
important protocol missing in our toolbox is a feasible 
distillation protocol for continuous variables. We are 
particularly interested in distillation protocols for entan- 
gled Gaussian states because these states can easily be 
generated in the laboratory. The protocols suggested so 
far involve rather complicated nonlinear transformations 
such as subtraction of a single photon M or quantum 
non-demolition measurement of total photon number in 
several modes Q . It would be of great help to have a dis- 
tillation protocol for continuous variables that could be 
implemented with linear optics and which would distill 
Gaussian entangled states. However, no such protocol is 
currently known and it is an open question whether such 
distillation protocol exists at all. 

In the present paper we attempt to shed some light 
on this issue making use of the formalism of Gaussian 
completely positive (CP) maps jlO[ O]. These maps rep- 
resent all transformations that can be carried out with 
the help of passive and active linear optical elements, ho- 
modyne detectors, and auxiliary optical modes prepared 



initially in Gaussian states. These transformations may 
be deterministic or probabilistic. In the latter case we 
accept or reject the output state in dependence on the 
output of a quantum measurement (with some Gaussian 
probability distribution). We shall consider an arbitrary 
bipartite probabilistic Gaussian operation which can be 
implemented with the help of local operations and clas- 
sical communication (LOCC). We shall prove that for 
input bipartite Gaussian states it is always possible to 
construct a deterministic LOCC Gaussian transforma- 
tion that yields the same output state (for a fixed input) 
as a given probabilistic LOCC Gaussian transformation. 
This implies that it is impossible to distill more entan- 
glement from a single copy of entangled Gaussian state 
by means of Gaussian operations. This should be con- 
trasted with distillation protocols for a single copy of 
two-qubit entangled state where the LOCC operations 
may in some cases allow to extract more entanglement 
p8fl . In particular, any pure entangled two-qubit state 
can be transformed with certain probability via LOCC 
operations onto maximally entangled Bell state. Further- 
more, our results imply that any hypothetical Gaussian 
distillation protocol optimized for given shared entangled 
Gaussian states would be a deterministic protocol and we 
shall find a generic structure of this optimum protocol. 
A version of this optimum protocol where Alice and Bob 
share two identical copies of Gaussian state with sym- 
metric covariance matrix has been considered by Eisert 
et al. who employed the log-negativity as the entan- 
glement measure and proved that it is impossible to dis- 
till entanglement via this protocol. These findings thus 
strongly support the conjecture that it is impossible to 
distill entangled Gaussian states via Gaussian operations 
and that some nonlincarity is necessary 

We shall extensively exploit the Jamiolkowski iso- 
morphism fl3ll between completely positive maps M. 
and positive-semidefinite operators (bipartite quantum 
states) x on tensor product of input and output Hilbert 
spaces TL®K,. In terms of the operator x > the rela- 
tion between input and output density matrices can be 
written as a partial trace over the input space, 



Pout = Tr in [xp£ ® H-out], 



(1) 



2 



where T stands for the transposition in some fixed basis 
and l ou t denotes an identity operator on output space. 
The operator x can be obtained from a maximally en- 
tangled state on H® 2 , |V> = £j =1 |i)iU) 2 , (d - drrnW) 
if the map A4 is applied to one part of this state, 



X = M®I[W){il>\] 



(2) 



Here X stands for the identity transformation. 

In continuous variable systems, we deal with infinite 
dimensional Hilbert spaces and the maximally entangled 
state \ip) becomes a tensor product of N in (unphysical) 
two-mode infinitely squeezed vacuum states, where _/V; n 
is the number of input modes. Gaussian completely pos- 
itive maps are defined as maps which transform Gaus- 
sian states into Gaussian states. Gaussian CP maps are 
thus isomorphic to bipartite Gaussian quantum states 
X- Now any Gaussian state \ is completely character- 
ized by the first and second moments: mean values of 
quadratures and a covariance matrix F. Define vector of 
quadratures r — (x\,pi, . . . ,xn,Pn) T where N is the to- 
tal number of input+output modes. The elements of ma- 
trix r are defined as Yij — (Ar^Ar.,-) + (Ar^-Ar*), where 
Ar^ = ri — (r<). Nonzero mean values of the quadra- 
tures of the Gaussian state x representing a Gaussian 
CP map indicate that this map involves certain displace- 
ments. However, these operations can be performed lo- 
cally and are therefore irrelevant for the entanglement 
properties and can be omitted. Thus we can assume that 
(fj) = and the CP map x IS lun y described by the co- 
variance matrix T. It is convenient to split the matrix T 
into input and output parts and write 



(3) 



where A stands for the covariance matrix of the "in- 
put" modes, B is the covariance matrix of "output" 
modes and C contains the input-output correlations. The 
input-output transformation rewritten in terms of the 
Wigner functions reads 

/•OO 

W out (r out ) = {2w) Ni ° / W x (r in , r out )W in (Rn n ) dr in , 



(4) 

where R = diag(l, —1, 1, —1, . . . , 1, —1) is a diagonal ma- 
trix which represents the transposition in phase space 
(xj — > Xj, pj — > —pj). It is convenient to deal with char- 
acteristic functions which are Fourier transforms of the 
Wigner functions, 



W(f) exp (ir ■ q) dr. 



(5) 



On expressing all Wigner functions in terms of the char- 
acteristic functions, we obtain from Eq. (0), 



C(q 



(27T) 



C x (<fin, qout)Cin(-Rqhi) dq u 



Assuming input Gaussian state with covariance matrix 
Tin, Ci n (g) = exp (— \q r T m q) , we find that the the out- 
put state is also Gaussian with covariance matrix given 
by 



Tout = B - C T (A + RTinR^C. 



(7) 



We now prove a very important feature of Gaussian 
CP maps. It holds that for every input Gaussian state 
and a probabilistic (trace-decreasing) LOCC Gaussian 
CP map there exists a deterministic (trace-preserving) 
LOCC Gaussian CP map which transforms the input 
state into the output state with the covariance matrix 
([?]). The explicit construction of the trace-preserving 
map is inspired by recent works on the possibility of stor- 
ing quantum dynamics in quantum states |l5). The 
basic strategy is to encode the transformation into a bi- 
partite state x which then serves as a quantum channel 
in the teleportation. In this way, the desired transforma- 
tion is carried out with certain probability depending on 
the dimension of the Hilbert space. 

The continuous-variable analogue of this scheme goes 
as follows. We prepare a Gaussian state x with covari- 
ance matrix T given by Eq. (^) and carry out a Bell 
measurement on the input state and the input modes 
of the state x- This measurement is performed sepa- 
rately for each corresponding pair of modes and consists 
of measuring the difference of x quadratures and sum of 
p quadratures by means of homodyne detectors. Let the 
vectors Xd and p^ contain the measurement outcomes for 
x and p quadratures, respectively. The (non-normalized) 
Wigner function of the output modes conditioned on the 
measurement outcome fd reads 



Wo U t(r ou t|r*i) = / W x (r in ,r out )W in (r)S(xi n - x ~ x d ) 

J —OO 

xS(Phi + p-pd) dx m dp in dx dp. 

Consider input Gaussian state. It turns out that the 
covariance matrix of the output state is given by Eq. (0) 
and does not depend on the measured quadratures Xd 
and pd- However, the output state is displaced by 



f CO nd = C T {A + RT in R)- 1 



(8) 



(6) 



If we know the input state and the transformation, then 
we can calculate f CO nd for given detected quadratures Xd, 
Pd and by means of suitable displacement transformation 
applied to output state we can always set the coherent 
signal in the output state to zero. In this way we obtain 
in a deterministic manner an output state which has the 
covariance matrix (0). This procedure is essentially the 
Braunstcin-Kimblc scheme for teleportation of continu- 
ous variables 0. We should note here that many Gaus- 
sian CP maps, in particular all trace-preserving maps, are 
represented by unphysical (infinitely squeezed) states x- 
However, we can approximate such unphysical state by a 
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FIG. 1: Setup for implementation of deterministic LOCC 
Gaussian CP map that is for certain input Gaussian state 
equivalent to a given probabilistic LOCC Gaussian CP map. 



physical finitely squeezed state with an arbitrarily high 
accuracy and thus also approximate the transformation 
(0) with arbitrarily high precision. 

Let us now turn our attention to the LOCC Gaus- 
sian CP maps. Obviously, every LOCC Gaussian map 
•Mlocc is isomorphic to a Gaussian state \ which is 
separable with respect to Alice and Bob. (Note that ac- 
cording to Eq. (H|) Alice and Bob can prepare the state 
X in their labs via LOCC operations.) A scheme for de- 
terministic implementation of any LOCC Gaussian CP 
map (for known input Gaussian state) is shown in Fig. 
1 for the simplest case when there is a single input and 
a single output mode on each side. By means of LOCC 
operations, Alice and Bob prepare the four-mode state x 
representing the CP map (Al and Bl are input modes 
and A2 and B2 are output modes). Alice mixes her part 
of the input state in mode A- lri with Al on balanced beam 
splitter and measures xai — XAm and PAin+PAi by means 
of two balanced homodyne detectors. Bob performs the 
same operations with his modes Bl and B m . Alice and 
Bob exchange the results of their measurements via clas- 
sical communication channel and apply appropriate local 
displacement transformations to the modes A2 and B2 
thereby producing the two-mode output state A ou t and 
Bout- Notice that the protocol works only for Gaussian 
states. If the input state is not Gaussian, than it may 
happen that some trace-decreasing LOCC Gaussian CP 
maps will yield outputs that cannot be obtained with any 
trace-preserving LOCC Gaussian CP map. 

A very important implication concerning distillation 
protocols is that we cannot distill more entanglement 
from a single copy of two-mode bipartite entangled Gaus- 
sian state by means of LOCC Gaussian operations. This 
follows from the fact that any probabilistic LOCC Gaus- 
sian operation can be replaced by a deterministic one 
which yields the same output (for a given input Gaus- 
sian state). However, any reasonable measure of en- 
tanglement must be non-increasing under deterministic 
LOCC operations. Here we note that this impossibility 
was noted in a paper by Parker et al. |f(|, who con- 
sidered distillation via continuous-variable entanglement 
swapping. While this procedure works for Schrodinger 



cat states, it fails for two-mode squeezed vacuum. Parker 
et al. also noticed that in the latter case, the probabilis- 
tic Gaussian transformation becomes deterministic. Our 
arguments show that this is a general feature of Gaussian 
states and Gaussian CP maps. 

Let us have a more detailed look at the structure of 
bipartite LOCC Gaussian CP maps. It was shown by 
Werner and Wolf jlTj that a bipartite Gaussian state is 
separable if and only if it can be transformed via lo- 
cal symplectic transformations into state with positive 
Glauber-Sudarshan representation, i.e. a state which is 
a convex mixture of coherent states and is not squeezed. 
We can thus write 

/oo 
P(a, (3)S A \a)(a\S A ® S B \p)(f3\S f B dadfJ, 
-OO 

(9) 

where P(a, (3) > is classical Gaussian probability distri- 
bution and | a), \/3) denote (multimode) coherent states of 
Alice's and Bob's modes. We have seen that it is impos- 
sible to distill a single copy of a Gaussian state by means 
of LOCC Gaussian transformations. What if Alice and 
Bob posses several copies? Assume that Alice and Bob 
apply the LOCC Gaussian map (^) to their states. This 
distillation map takes all copies as an input and yields a 
single copy of two-mode state shared by Alice and Bob as 
the output. The output state is a mixture of states with 
identical covariance matrices and varying displacements. 
In terms of Wigner functions, we can write 

/oo 
P(a,f3)W(r-r d (a,f3)) da dp, (10) 
-oo 

where rd(a,(3) is the displacement and P(a,f3) > 0. 
However, all states with Wigner functions W(r — fj) ex- 
hibit the same entanglement, because entanglement de- 
pends only on the covariance matrix and not on the dis- 
placement. Hence it is always optimum to choose an 
LOCC CP map which is represented by a pure Gaus- 
sian state. An analogous situation arises in distillation 
of a single pair of entangled qubits where it is optimum 
to apply local filtering operation |Q (this is a trace- 
decreasing CP map whose Kraus decomposition contains 
only one term and the CP map is thus represented by a 
pure state). 

Since we assume that we know the state that we want 
to distill, we can transform any LOCC trace-decreasing 
map onto trace-preserving map hence the optimum pro- 
tocol will be deterministic. Consider now the simplest 
case when Alice and Bob share two pairs of entangled 
Gaussian states. The transformation on Alice's side is 
represented by a pure three-mode Gaussian state, which 
can be obtained from three-mode vacuum via some three- 
mode symplectic transformation. This three-mode state 
splits into two input modes and one output mode. By 
means of "local" symplectic transformations on input and 
output modes, we can transform the covariance matrix 
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FIG. 2: Optimal distillation protocol for two copies of two- 
mode entangled Gaussian states. 



of this state into the form 
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(11) 



The reduced density matrix of each input and output 
mode is density matrix of thermal state and there are 
no correlations between the two input modes. Since the 
whole three-mode state is pure, the density matrix of the 
two input states has the same eigenvalues as the reduced 
density matrix of the output state. This is possible only 
if one of the input modes is in pure vacuum state. This 
leads to a = c, b = 1, e\ = e<i = e$ = 0. One of the 
input modes is effectively decoupled and we end up with 
pure two-mode squeezed vacuum state. We thus have 
a very appealing and intuitive picture: any pure three- 
mode state can be prepared if we start from (suitably 
chosen) two-mode squeezed vacuum state and a vacuum 
state and apply single-mode symplectic transformation 
to the output mode and two-mode symplectic transfor- 
mation Si to two input modes. We may thus write the 
transformation on Alice's side in the form 



XA = Si(XA0®\0)(0\)$, 



(12) 



where xao is an operator on Hilbert space of two modes 
(one input and one output). To see what are the impli- 
cations, we insert this expression into formula ([!]) and for 
the moment do not consider Bob's states. We get 



Pout 



= Ti'i 
= Tr ; 



Axao® \0)(0\Slpf n Si® lout 
1 [ X Ao®|0)(0|(5f Pin 5*) T 



(13) 



From this formula we can see that the transformation 
reduces to the following three steps: (i) apply symplec- 
tic transformation Sf to the input two-mode state, (ii) 
project the second mode onto vacuum state, (iii) Apply 
a CP map xao to the first mode. The transformation 
on Bob's side has the same structure. This protocol can 
be further simplified. After Alice and Bob project one 
of modes onto vacuum state, they possess only a single 
mode each. The application of the local maps xao (xbo) 
cannot increase entanglement, because we have shown 



that it is impossible to distill a single copy of two-mode 
entangled Gaussian state by means of LOCC Gaussian 
operations. Thus we need not consider the transforma- 
tions xao (xbo) and the resulting optimum simplified dis- 
tillation scheme is shown in Fig. 2. Alice and Bob share 
two pairs of Gaussian states. Alice's and Bob's modes are 
labeled as Al, A2 and Bl, B2, respectively. Both Alice 
and Bob locally apply some two-mode symplectic trans- 
formations Sa and Sb to their modes. Subsequently they 
both feed the modes A2 and B2 to eight-port homodyne 
detectors (EHD) thereby projecting them into coherent 
states \a) and |/3). Finally, they exchange the results of 
their measurements and displace appropriately the out- 
put states. This scheme is deterministic and represents 
the optimum Gaussian distillation protocol for Gaussian 
states. Eisert et al. proved that it is impossible to dis- 
till entanglement from two identical copies of two-mode 
Gaussian state with symmetrical covariance matrix via 
this protocol 12 . All these results strongly support the 



conjecture that it is impossible to distill entangled Gaus- 
sian states with Gaussian operations. 

Note added: After this work was completed, I learned 
that Giedke and Cirac (rJ have also investigated the 
properties of trace-decreasing Gaussian CP maps and 
they independently obtained similar results. Moreover, 
they proved that the distillation of Gaussian states with 
Gaussian operations is impossible for an arbitrary num- 
ber of modes per site. 
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